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I. INTRODUCTION 


A. BACKGROUND 


In today's world, the ocean scientist and engineer need 
increased accuracy for marine positioning to accomplish 
their work. The establishment of a system for the precise 
determination of positions on the deep seafloor is one of 
the most fundamental challenges of working in the open sea. 

A seafloor benchmark positioning system is proposed 
which can determine accurately the location of- objects on 
the sea floor (Saxena, 1974). Such precisely located bench- 
marks can later be used in turn to delineate offshore prop= 
erty lines and national boundaries. They are also useful in 
solving problems in seafloor engineering, in plate tectonic 
studies and in connection with bore-hole reentry associated 
With off-shore oil recovery. 

The objective of this study is to demonstrate the feasi- 


bility of establishing such seafloor benchmarks. 


B. DESCRIPTION 


Early experiments associated witn seafloor acoustic 
transponaer arrays were carried out by Hart (1967), Haehnle 
(1907), Fubara and Mourad (1972), and McKeown (1975), who 
used a baseline crossing method to soive for the distances 
between transponders. "Baseline" in this thesis means the 
line between the projections of two bottom-mounted acoustic 
transponders onto the sea surface. 

Ihe baseline accuracies obtained by Hart (1967), McKeown 
(1975), and Fubara and Mourad (1972) were +15.7 m for base- 
lines of 5509 m, +3.1 to +3.2 m over baselines of 6373 to 


7219 m and +15.5 m over baselines of 9364 m, respectively. 
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The first two of the above authors did not use a least 
squares method for data analysis. 

There are three parts to the determination of the abso- 
lute position of an ocean-bottom transponder array using 
satellites. They are (Figure 1.1): 

1. То determine the geodetic position of tne ship's 
receiving antenna from the Global Positioning System 
{GPS) ; 

2. To convert the geodetic position of the antenna to 
the ship- mounted acoustic transducer; anā 

3. To determine the geodetic position of tne acoustic 
transponder array on the sea floor using a ship- 
mounted acoustic transducer. | 

This thesis is concerned with the relative position oí a 


bottom-mounted transponder array with regard to tne ship. 


C. DESCRIPTION OF THE SYSTEM 


Acoustic transponders manufactured by Oceano Instruments 
were used for this project. They were deployed by the R/V 
Acania. Figure 1.2 shows the configuration or the Oceano 
system used. Their principal features are: 
1. АМ-121 : (Acoustic Module) 
The hydrophone was fixed amidships to the 
port Side at a depth 2.95 m. 
- Transmission irequencies : 8 to 16 kHz. 
- Reception frequencies : 8 to 16 kHz. 

2. 851-121 : (Recoverable Transponder) 
- Frequencies : 8 to 16 kHz. 
- Delay time : 15.0£0.1 ms. 

3. 11-201 : (The Accoustic Telecommaud Moauie) 
has the following functions: 
- Code acoustic signal transmitted by either tne 

acoustic module or UQC transducer. 


- Monitors range to desired transponder. 
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The Configuration of Bottom Transponder Array 


Shown with Ship and GPS Satellite. 


Figure 1.1 
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- Monitors reception of code by transponder. 
Ц. ІМ-100 : (баға interface Module) 
- Accepts up to 16 data inputs. 
- Max. converted data format : 7 digits without 
Sign or 6 digits plus sign. 
- Computer interrupt request possible through 
any input. | 
- Each channel may be read separately. 
5. | RM-201 : (Rangenmeter) 


Transponder turnaround delay offset : 0 to 99.99 ms. 


- Mean velocity of sound input : 0 to 9999.9 units per 


second. 


Standard reception filters : 8 to 16 kHz (0.5 kHz step) 
Bandwidta : 200 Hz. 


Critical to the experiment are the following: 


1. The transducer (а part of the AM-121) was 2.95 m 
below the surface, and the transmission and reception 
rrequencies were 8 to 16 kHz. 

2. The system delay on the RM-201 was set at 15 miilise- 
conds,  i.e., the time iapse between reception and 
transmission of a signal. 

3. Velocity on the RM-201 was set to 1480 m/s. 

4.  Interrogation period was set at 10 5. 

5. Interrogation frequency (kHz) was set 15 kHz. 

The numerical constants of (3) and (4) were used for the 
survey. Details of the transponders are given py Oceano 
Instruments (1984). 
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A. SITE SELECTION 


The main criterion for tne site óf the experiment was 
that it be relatively flat over an area approximately 2 km x 
2 kn. Since there were no existing cnarts detaiied enough, 
a reconnaissance survey was perfomed by the R/V ACANIA on 10 
and 11 April 1985. Some 10 north-south sounding lines and 
l-cross iine were run with 1-nml spacings and 5-min position 
fixes. 

The area chosen for the experiment (Figure 2.1) is 
approximately 24 nmi west of Pt. Lobos, California, in water 
approximately 1500 to 2000 m deep. It is the area from 
29.5" to 36931'N and 122916' to 122917.8'H. Mini-Ranger 
Stations (Figure 2.2) were set up at Point Sur 
feo 181 20.279" Н X 1219 53°56. 179"W) and Carmel 
(36933'49.176" N X 121953'48.358'"W) for additional ship 


position control. 


B. | IMPLANTATION PROCEDURE 


The transponders were deployed at pre-selected sites 
using a buoy-first, anchor-last, technique. The anchors 
used were railroad wheels which weigh about 290 kg in water. 
The mooring cables between the anchors and the transponders 
and between the transponders and the buoys are each 8-m 
long. The transponders weigh about 46 lbs in water, and 
each buoy has 55 lbs positive floatation. Three 17-in 
Benthos glass sphere floatation buoys were attached to each 
transponder by meats of 5/32-in stainless cables 8-m long. 


Altnough 2 spheres would have provided sufficent buoyancy 
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for recovery, a third was used to provide additional tension 
(about 50 kg total) to reduce the motion of the transpon- 


ders. Тһе mooring system is shown in Figure 2.3. 
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Figure 2.3 Schematic of Bottom Transponder Mooring. 
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III. DETERMINATION OF HARMONIC MEAN SOUND VELOCITY 


A. INTRODUCTION 


Historically, echo soundings have been made by assuming 
a constant, approximately average, value for the velocity of 
sound throughout the water column, usually 
4800 ft/s (21463 п/5) or 1500 m/s, and then correcting for 
variations of this assumed value for the actual water 
column. Echo sounders are time-measuring devices. The 
sounder's acoustic transducer emits a sound pulse which is 
reflected upwards by the bottom-mounted acoustic transponder 
and received back at the surface transducer. The time of 
travel of the sound pulse is divided by two, and this value 
is multiplied by the assumed harmonic mean speed of sound in 
sea water, thus, giving the distance according to tne 
expression D - V t . The echo sounder makes this transfor- 
mation electronically or mechanically within the device 
itself and displays the aistance in the water between trans- 
ducer and transponders. This distance is not equal to the 
true distance, since the assumed harmonic mean sound speed 
generally does not equal the true mean speed for a partic- 
uiar depth. To determine the true harmonic mean sound speed 
one must know the sound speed throughout the water coiumn 
and apply it to solve for the harmonic mean. For a distance 
measurement to be accurate, precise measurements of time and 
the harmonic mean speed of sound in tne water column are 


mandatory. 


B. THE VELOCITY OF SOUND IN THE OCEAN 


The velocity of sound in sea water is a function of 


temperature, pressure, and salinity: V = V(T,S,P). Because 


20 


tne sound velocity is not constant with position, Sound rays 
are refracted according to Snell's law. 

We have followed the standard practice ol the U.S. 
National Ocean Service (NOS) to use Wilson's (1960) equation 
for the speed of sound as a function of salinity, tempera- 


ture and pressure (Umbach, 1976): 


3.1 
V21449. 14* V, * V5 * Vs * V тр m 
where 
V,24.5721T-4.4532X107?17-2.6045X10 ^T^ 
(3.2) 
+7.9851x10 *T* 
220 -5.2 
Vp=1.60272X10""P+1.0268X10""P n. 
xg: -12.4 
+3.5216х1079р3-3.3603х107 ‘р 
Vy=1.39799(5-35)+1.69202X10""(5-35)* (3-4) 
Vomp=(S-35)(-1.1244x10"%T+7,7711x107 T 
+7.7016X10"%p-1.2943X10'P*+3.1580x10"*pT 
+1.5790Х10 ?PT?j«P(-1.8607X10 ^T 
+7 4812Х10 91?«4.5283x10 ^91) | (3.5) 


7 


+P2(-2.5294x10""T+1.8563x10""T*) 


+P%-1.9646x10"""7). 
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In these equations the absolute pressure, P, is expressed in 
kg/cn?, temperature, T, is in ?C, salinity, S, is in g/kg, 
and sound velocity, V, is in m/s. For these relations for V 
the standard deviation from the mean is 0.30 m/s for all 
data obtained in the ranges -4 < T< 30 9C, 1 < P < 1000 
kg/cn*, and 0 £95 < 9 J/kg ARIT Ssona, 1960, PAION 


С. HARMONIC MEAN SOUND VELOCITY 


The slant range from a ship to a transponder is deter- 
mined by measuring the transit time, t, of an acoustic puise 
and converting it to distance using an appropriate value гог 
sound speed. If the mean sound velocity, V, in the water 
column is known, the distance, L, from transducer to tran- 


sponder can be computed by: 
һ - (ЕН (3.6) 


For our small project area the distances between transpon- 
ders are less than 2 nmi, and it is assumed that there is no 
hocizontal sound speed variation, and oniy corrections for 
vertical sound speed variation are considered. For a 
one-way travel time , T, through a water column of depth, 2, 
the mean sound velocity is: 


У = ИТ (3-7) 


The meaning of "mean sound velocity" is neither "mean 
velocity trom the surface to the stated depth" nor the "nean 
value for the velocity of sound through the vertical water 
column". According to Equation (3.7) the mean velocity is 
obtained from tne true depth divided by the one way travel 


time. 
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Figure 3.1 A Layered Ocean Where Y; cat the Sound 


= ү І 
Speed in the ith Layer and Z IS © 


2) 5 
Depth. 
Consider now a layered ocean having a different sound 
Speed in eaca layer. Figure 3.1 shows a series of finite 
layers of thickness А2., each with an associated sound 
speed, Vi - The time interval, AT., for the sound wave to 


pass vertically through the ith layer is: 


AT: - AZ; / NV (3. 8) 


Sunning tne time intervals for all n layers in a vertical 


water column from surface transducer to bottom transponder , 
yt n 

= > AT. = E (Az; /V ) (3. 9) 
=, p) 


For a continuous function, V - V(z), this becoues 
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4 
T = | dz/V(Z) (3. 10) 


о 


Equation (3.7) may now be written : 


7 = 24 (32 485 


: dz 
о УЕ) 


This is the integral form of the harmonic mean. Using 





Simpson's Rule, Equation (3.10) can be rewritten as : 


т = (ЛУЧ РАТНЕ пече +17V_) 2/3 (3. 12] 
where d - D/n (magnitude ot depth increments) 
р = pottom depth at transponder 
n = even number of depth increments between 
transducer and transponder 
ү; = measured sound speed at depth j. 
Since the harmonic mean sound velocity is obtained by 
Equation (3.11), the slant range from transducer to acoustic 
transponder can be obtained using Equations (3.10) and 


(a. БРИЈЕ 
AA TU. 


СТ) casts to 1400 m were made on 18 and 20 Мау 1985 іг 
tae project area just prior to and during our main experi- 
ment. Sound velocities and harmonic means based on the 
sound velocity profiies derived from the CTD casts are 
plotted in Figures 3.2 апа 3.3. Although the profiles 
Showh in these figures are for different days and different 
places within the project area, they Show very little varia- 
tion in harmonic Mean sound velocity, and thus support our 
previous assumption tnat the sea water is horizontally homo- 


geneous to an acceptable degree within the study area. 
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HARMONIC MEAN COMPARISON 
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Figure 3.4 Harmonic Mean Comparison. 
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Actually, the difference of these two harmonic mean veloci- 
ties is approximately 0.07 m/s near the bottom as is shown 
in Figure 3.4. The mean of these two narmonic mean sound 

velocity profiles is used in the range calculations of the 


zollowing chapter. 


D. ACOUSTIC WAVE REFRACTION 


To simplify our acoustic analysis it is desirable to 
assume that in the study area sound follows straight paths 
from the surface to the ocean bottom. То test this approxi- 
mation the differences in path lengths from the surface to 
the bcttom were calculated for varies depression angles at 
the surface. 

Profiles such as those of Figures 3.2 and 3.3 are 
usually simplified for analysis by separation into an appro- 
priate number of segments each having an approximately 
constant gradient (Figure 3.5). If the velocity of sound 
changes linearly with depth, sound rays can be shown to nave 
a constant radius of curvature (Urick, 1983, p. 124). 
Horizontal distance, d, depth, z, radius of curvature, R, 
segments of arc, S, and depression angle, A, arce shown sche- 
matically in Figure 3.6. The path of a ray when the speed 
or sound varies witn depth can be calculated by application 


of Snell's law: 
cosA, / С, = COSA, IE EG совА./ С, (3.13) 


where C;= sound speed 
A: = depression angle 
For an initial depression angle oí 30% the profile of 
Figure 3.5 gives the results shown in Table I. 
Ihe raaius of curvature can be obtained in the following 
Manner (Kinsler, 1982, pp399-402). Let 
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| 
| 
| 
| 
| | Sound | Depth | Depression 
| Points | Velocity | (mn): | angle (deg.) 
| — енә но еә ШР ор НР Әәә а р «ә «ШР «ыы саш елә «әнә «ШР «иш «ӘР ә — ———— —— ——.———— ——— и oe es es о ша а со ар св = > 
| 1. | 1493.0 | 0. l 3029 
| zs | 1483.5 | 176.4 | SJØL GE: 
| За | 1484.04 | 266.6 | БВ 
4. | 1481.2 | 526.4 | 507 78 
4 Da | 1482.46 | 1097.2 | 30569 
| б. 1 1484.34 | 1400.0 | 30257 
E - B 

ape = ( E 571727 zn 2. ) (3. 14) 
and 

B; = - (С; / 9. } COSA. (3. 15) 
where 


The change in range аа. and depth 47; аге (Kinsler, 1982) : 


Az. = ~R: ( COSA, 


TABLE I 
Depression Angle at Different Layers 


R; = Radius of curvature in the ith depth increment 
C; = Sound speed in the ith depth region 
A; = Depression angle 


g. = Speed gradient 


да. т FRH ( sinA; - S1IDA¿,,) 


С L ] 
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(3.16) 


= COSA, ) (35807) 
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so the chord can be obtained by two components of change in 
range and depth; that is, 


D. = ( ad2 + А22 )1/2 (3. 18! 
L L С 
The angle, 9: , 15 
0, = 251п-1 ( D; /2R; ) (3. 19) 
so the arc, S, gs. 
5. = К. Ө. (3.20% 


where 9: 


The range and depth increments give tne straight line 


is the angie in radians. 


ШЕИ ОСИ ло куе е (32-21 
The ray patis can be obtained by summing arcs: 
S = а (3.22 


The difference due to refraction between a sound ray's 


actual path and a straight line is thus: 
Es = i (3:23 


where Е - the error between two lines, 


the distance along the ray path, and 


it 


the straight-line distance between the 

Ship-mounted transducer and bottom transponder. 
Ray-trace computations were made for various initial 

depression angles (at the surface) down to a depth of 1400 m 


(Tables 11 and [111);Since E<1 cu for | - пл ас 


D 


As 


Tne Error 


TABLE ІІ 


The Result of Ray Trace Computation 
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For Initial Depression Angle = 30 Degrees 


= oe ee ee ee ee aa ee ee ee ee ee ee TO TA === 


0.0096 meters 


Depth Gradient Radius porn: Depth Arc 
Increment g R Dist. (å ) 2. 5 

w o ooo In)... - в ml. ml А 
| illa 0.0538 320 м. 36 301.72 176.4 349. 5023 
е 0.0060 287966.49 152.47 2052 17721553 
| 3. 0.0109 157706.64 437.85 259.8 509.1278 
| 4. 0.0022 750994782 960.03 57028 1116-8993 
| 5; 9006 2277609187 511.35 3012: 8 594.2820 
| Tot: 2363.42 1400.0 2746.9658 
| The Distance of Straight Line : 2146. 9562 
| 

| 
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depression angle and a 1400-m water depth, corresponding io 


a horizontal distance of 2363 m, 


which covers the project 


area, and since the depression angles encountered during the 


project were 30% or greater, 


1 cm represents an upper limit 


to the error resulting from our assumption of linear ray 


paths from the surface to the bottom, and a correction for 


refraction is unnecessary. 
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TABLE III 
The Error Comparation at Different Angles 


Depth : 1400 meters 


=> => «=> E A A A A ee A қы» «ы» «қы» «ы» кы» «ғы» оғы» O A «> «ap A A CA ЕР ЕРЕ ЕР о TO «ғыр «кы» «Жы» «кы» «кы» «кы» GOE A «ED НР НР A: CURED «нә» 


Depression  Horizantal Ray Straight Error 
Angle Distance Trace Line 

(degree) (n) (m) (n) (са) 
0 134529223 13508.3850 13502.0064 637.86 
10 6637.03 6783.2984 679320595 21789 
20 3647.10 3906.6140 3906.5806 3:38 
30 2363.42 2746. 9658 2746.9562 Gn Fo 
60 801. 28 1613450859 1613.0853 9.06 

90 0.0 1400.0 1400.0 0. 0 


34 


L 


IV. RELATIVE POSITION DETERMINATION OF IHE TRANSPONDER ARRAY 


A. INTRODUCTION 


Two methods to determine relative positions are 
discussed in this chapter. They are basic baseline crossing 
and the modified baseline crossing method; the former is 
described by Hart (1967) and Haehnle (1967), the latter is 
described by Fubara and Mourad (1972). The main difference 
between the methods is that Hart and Haehnle require accu- 
rate transponder depths, but Fubara and Mourad- do not. 

The cloverleaf method used to determine transponder 
depths 1s described by Haehnle (1967) is introduced in tae 


next section. 


В. CLOVERLEAF METHOD FOR TRANSPONDER DEPTH DETERMINATION 


An initial approximation to the depths of bottom-aounted 
acoustic transponders can be obtained from a knowledge of 
the water depth at the point where the mooring is released 
from the ship. However, the unit may nove laterally while 
descending, or bottom topography may have considerable 
reliet not evident on the echo sounder record. Hence, a 
more accurate metnod of depth determination is required. 
Here the depth is defined as the distance between the huli- 
nounted trarsducer and the bottom-anchored transponder waen 
the vessel is positioned exactly over the transponder. 

The ship need not be directly above the transponder buat 
Should be in close proximity to it. Although it is very 
difficult to cross over the transponder exactly with tne 
Ship, the depth may be approximated as tne minimum slant 
range when the transponder is crossed over many times using 


a "cloverleaf" method (Figure 4.1). Thus, if a transponder 


B 


=) 


transponder (drop point) 


| A location of the 


a. The first pass through the transponder location. - 


A “ 
| \ 
GE \ 
Y 13 Drop Point 
,Fourth CPA 
Pirst CPA | 
9.909. со а 

> — TT т ш + 

= AN “Third CPA 
Secony CPA 


b. Maneuver for obtaining four CPA's. 


OEA, DD COHEN 


Figure 4.1 Determination of Transponder Depth. 
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Figure 4.2 Exarple of Depth Difference. 


is at a depta of 1500 m and tae Ship 1S displaced 50 meters 


horizontally from it, tne depth determined by the cloverieaz 


method would be in error by only 0.33 ueters (Figure 4.2). 


The procedure used for running a cloverleaf maneuver is 


as follows (Haennle, 1957): 


КЕ 


2. 


Each time a transponder is dropped from the ship the 
Ship's position is noted. 

The ship proceeds toward the transponder drop point 
and passes near it; А closest point of approacn 

(CPA) is then determined and recorded. 

МАО пен turns. crossinj the first track perpen- 
dicularly and passing near the first CPA. Once 
again, the CPA of this run is recorded. 

Tne above maneuver is repeated until at least four 
such CPA's nave been obtained, each time attempting 
to pass through the previous CPA perpendicular to the 


track. 


SÅ 
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Pigure 4.3 Sample of Depth Determination. 


5. The depth of the transponder 1s determined by plot- 
ting the slant ranges versus time for each of the 
гоог CPA runs. The minimum point of tais curve is 
considered as the minimum slant range for that 
particular transponder (Figure 4.3). 

The R/V ACANIA ran Siwliar cloverieaf patterns оп 21 Мау 
1985. She crossed near each transponder many times as 
closely as possible. The minimum distances and corrected 


depths are listed in Table IV. 
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TABLE IV 
Transponder Depth by Cloverleaf Maneuver 


(From Oceano [ 1985, p.-12]) 
Corrected 
Transponder Dept h 
(n) 

1366.9 
1394 . 3 
ЗӨ 55 
1406.7 
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C. BASELINE CROSSING METHODS 


Baseline length can be determined utilizing tne baseline 
crossing method described by Haehnle (1967) and Hart (1967), 
which involves repeatedly transiting the baseline while 
simultaneously measuring ranges to two transponders. The 
Minimum sum of the two horizontal ranges as determined fron 
measured slant ranges and known transponder depths is taxen 
to be the baseline length. 

The basic baseline crossing method reguires a kaowledge 
ОЁ: 

1. The ship's course and speed, 

2. The accoustic slant range (or the two-way travel 
times from the ship to two adjoining ocean- 
bottom transponders, and 


3. The harmonic mean sound velocity. 


Ss 





Sample of the Baseline Crossinge 


Figure 4.4 
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The ship, travelling at a constant speed and heading, 
crosses the lines joining the transponders as shown in 
Figure 4.4. The acoustic slant ranges are recorded as func- 
tions of time during the crossings. Тпе sum of two hori- 
zontal ranges (the distance T1-T2 in Figure 4.4) is minimum 
when the ship is in the vertical piane containing the two 
transponders. Similarly, for a four-transponder array, the 
distances Т1-Т2, 12-13, Т3-Т4, Т1-14, T1I-T3 and T2-T4 are 
determined. The minimum sums can be obtained by least- 
Squares fitting the horizontai ranges against time and 
obtaining the vertex of the parabola as shown in Figure 4.5. 
During the experiment four to six crossings were made for 
each baseline between transponders. The results are 
discussed below in Sections 1 and 2. 

Although the baseline crosSing method is relatively 
Simple with respect to data reduction, it is operationally 
time consuming. During the experiment it took one nour and 
toirty minutes for each baseline crossing. Thus the eight 
crossings (4 for north-south crossing, 4 for east-west 


crossing) took a total of 12 hours. 


1. Basic Baseline Crossing Model 





Figure 4.4 shows a Ship track crossing a baseline. 
The ship's transducer was set to interrogate the bottom 
acoustic transponders every 10 seconds for this experiment. 
Ranges were obtained from each transponder. It is easily 
seen (Figure 4.4) that, as the ship approaches the paseline 
TiT2, the total horzontal range at the first ship point (0) 
is greater than at the second ship point (1), and so on, 
until a minimum sum is reached at the instant of baseline 
crossiag; from there on the total increases for each pcint 
on the track. Since the ship maintains a constant course 
and speed, the mathematical function which describes this 
relationship is a parabola (Figure 4.5) of the form (Hart, 
28967): 
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Figure 4.5 Sample of Range-Time Parabola Showing 
baseline Crossing Tine. 


R=AT2 + BT+C (8.1) 


wnere 
A, B, C are constants to be determined, 
15 tine [s], and 
is total horizontal range between transponders 


pmj. 


To úeteraine the unknown variables A, B and C in 
Equation (4.1) a least squares adjustment is performed. At 
least four sets of data К = R(1) are needed to solve 
Equation (4.1). Since observed parameters contain randon 


error, Equation (58.1) can be rewritten as: 
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AO ABACO = * V (4. 2) 


where 


V is the residual of the observation. 


In matrix notation this corresponds to: 
2 2 
Vi -T, , ЕЛІ! A | К, 
5 | 
V2 -Tz pom Т; / = А Ка 


ЖҰ” A | (4.3) 
| | 


Been Hd 


which may be written as: 


He cT X, =R MEE 
MK! | x3 3X1 NI 


where LA are the residuals of the observations, 
I. are the numerical coefficients of the 
unknown variable, 
X; are the unknown variables(A, B, C), and 
R: are constant terms of total horizontal range. 


The standard deviation or the mean of sound velocity 
is +0.3 m/s (Wilson, 1960). So, the standard deviation of 
the distance is the following, where time is assumed error- 


less : 


Each observation is considered to be independent and 
uncorrelated. Since a measurement of high precision has a 
small variance, and one of low precision has a large vari- 
ance, the higher the weight the higher is the precision and 


vice versa. Accordingly, the weight W of a single obsera- 
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tion is defined as a quantity that is inversely proportional 
to the variances of the observations (62) (Mikhail, 1981, 
P2660). 


н = бг у 02 (ч.6) 
where 02 - variance of unit weight 
For the adjustment of Equation (5.4), the least squares 


criterion requires minimization of the weighted function G 
(Mikhail, 1981, pp.69-73): 


_ № 
$ = W VZ +W V2 + eee +R V2 = SH. V2 - (8.7) 


Loi SOE 


which in matrix form is written as 


ON Vu V (4.0) 
where 
woe 
OW 
и 


Wan 
Now Equation (4.4) can be rearranged as: 


V=R-TX (4.9) 


Substituting for V into Equation (4.8), get 


Q2 
| 


(R-TX)* W (R- TX) = (RS - X 1^) ui (8- TX) (4. 10) 
SO, 


Cae 
B* WR-X T Wh-B RTX*X T WTX (4. 11) 


42 


44 


Since G is a scalar, the right hand side of Eguation (4.10) 
is also a scalar. Furthermore, the transpose of a scalar is 


equal to the scalar: 


xT wR = (KUTU NB) = Вент (172) 


and 


чм $ (4. 13) 


Thus, Equation (4.10) can be written as: 


б = E* WR E 28* WTX * x* (T^ WT) X (4. 15) 


In Equation (4.14) all matrices are constants, except X, the 
matrix of unknowns. For G to pe a minimum, its partial 


derivative with respect to X must be equal to zero: 


2G > _овбут + 2х (ТнТ) = 0 (4. 15) 
ох 
and, 
E 6 
вн = W Е 
эх, хп дз) я Bon HAN : (Geo) 


giving the solution for X: 
Х = (rt W T )-1 TE yr (4. 17) 


Now, the minimum distance occurs when dR/dT = 0 


which is given below: 


SI O Å (4. 13) 
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Thus, to get the minimum distance of the sum, one may 
substitute the time from Equation (54.18) into Equation 
(121) - 


R (min.) = B2/4A - B2/2A ж С (4.19) 


«пеге A, B, C are given by Eguation (4.17). 


The variances of A, B and C (Mikhail, 1981, р258) 


are: 
t | 
92 - 02 (Т HT) -1 (4. 20) 
The standara deviation of T in Equation (4.18) 16: 


Ó- = [ (CG, B/2A2)2 + (6,722)? ]17? (4. 21) 
The standard deviation of the minimum distance R in Equation 
(4. 19) 18%: 


б, A Ip, yO A de (4.22 


2. Data Processing 


During the experiment at sea, both north-soutn anid 
east-west baseline crossings were run (Figures 4.6 and 4.7) 
using data collected by Oceano Instruments (1985a). 

For each baseline crossing, 11 consecutive interrogation 
points were fit to a parabola by least sguares. The reason 
for using 11 points is that this number is sufficient to 
define the parabola near its minimum. The minimum value of the 
parabola can be considered to define the baseline length. 


The parabolas shown іп the Figures 4.8, 4.9, 4.10, 4.11, 


46 


аараан PP D — ш. хш. шы Aa EE A EEE WA PP mm en En — al 


‚ТЗ 


ieee SS a SS Gs Se ee mue ro ERA is E A | 


1905 


T1 





2 


ГТ МЕНЕИ ——uuoUoÜÓ | 


equ eer que ec quo mna n qui aai ya - —p- -ң 
Бағы a EE lo жы сам сес ы Сары өле 


Figure 4.6 North-South Baseline Crossing. 


5.12 and 4.13 (at the end of this section) are tor single 
5345 between 11-12, 72-173, T3-I4, 71-74, Т1-Т3 апа 
12-75. Ine results are tabulated in detail in Table V, Which 
Gives the baseline lengths and their standard deviations ror 
each crossing. To get the mean of these lengths, one can 
use the foliowing equation: 

Em 7 (Eltk2* еее «Ел) / n = f(R) (4.23) 

Since the standard deviation of the mean length «cin 

Eguūuation (4.23) is: 


ба =0 (61 а /ак1)2+ (02 ӘҒ/Әк2)2% eee 11/2 (4.24) 


ме пауе 
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Figure 4.7 East-West Baseline Crossing. 
О =( (61/п) 2+(82/1) 2 + •• + (Оп/п) 2 J1/2 (4.25) 


Equations 4.23 and 4.25 are used with the data of 
Table V to compute tne mean of baSeline lengths and their 
standard deviations (Table VI). 

Lyman Burke oz Ocaano Instruments and L. Spielvogel 
or Seaco, Inc., have also computed the baseline lengtns. 
They use a different method and techrijue using most of the 
data collected during tne experiment. Their results and the 
differences between them are tabulated ir Table ViI- Due to 
our more limited data set and uncertainty in the depth for 
the baseline crossing method, these results differ fron 


tuose obtained with the crossing method by about 2 m. 


4 8 


PRA AAA RA AE AAA | 


TABLE V 


Baseline Lengths Determined by Baseline Crossings (m) 


Crossing | ПЕЕ T2 | T2 m.s | Tg mq 
| 1793.32:2.43 | 1993.09+40.381 |] 1972.39+2.50 
2a ГО 1991.92*2.73 | 1969.5042.02 
| | 
| | 


d 1793.6 7+2. 27 1992. 7542.38 | 1970.29+2.19 
1791.77+2.22 1994. 33+2.09 | 1969.97+2.23 

Crossing | 11 - T4 | T1 = T3 | 22 - TU 
12 1860.8 1+3. 25 | 2789.62+2. 20 258846642203 
De 1858.43+1.81 | 2789.77+1. 62 2589.70+1.94 


| | 

| | 
3. ЕСЕ ка „= | 2 09. 6542. 11 | 258748841. 14 
| 1858.48+2.08 | 2789.73+1.63 | 2588.17+2.21 

| | | 

| 


„у үштү тш „үк „ее үш кти ши ii сас)... T, fe RR e SA E 
e 
== ји и би Ба ду Mile cen. aman 


2. 18 58.6 5+2. 04 217 309722 t 2220 2799019. 7 2: 200 
6. | | 2580808282: 2.20 
i — == 1 





3. Modified Baseline Crossing Method 


Mourad and Fubara (1972) developed this method. 
Taeir method determines the relative positions of ocean- 
bottom transponders and the depth of each trarsponder. 
There is no necessity to accurately measure each depth 
directly; this is the main difference from the previous 
work. The physical principle of this method is shown in 
Figure 4.14 5 15 the baseline length and L1 and L2 are a 
pair cr coplanar ranges from the ship to transponders 1 and 
2 whose depths are Z1 and Z2, respectively. The mathmatical 


expression for this configuration is: 


49 





A A MS E A ыы: 


TABLE VI 
The Mean of The Baseline Lengths 


| Mean Value | Standard Deviation 


| ( m) | (а) 


AA a AA GE ME AA aa aa aa AA MA —— a Se O > AO осташе << AA WA Aa mug) сс => ===> == НЕЕ aa SE AAA O O ee 


¿A e a т AA aa O ai cc жалы (ығы 


71-72 | 1792. 84 | 1. 14 
T2 - T3 | 1992. 42 | 0.93 
T3 > T4 | 1970. 54 | Va 
T1- тұ | 1859. 00 | 1.04 
11-13 | 2789. 76 | 0.89 
T2 "PI 2588.74 | 0.82 
$ = {( 112 - 212 )1/2 + ( E22 - 722 у (4. 26) 


To determine the three unknowns S, Z1 and 22 by 
least squares method, at least 4 pairs of L1 and L2 are 


needed. Now, Equation (4.26) can be rewritten as foilows: 


-P ="S= (hile te yt/2 - ( 122 - 722 )1/2 (4.27) 
Assume 

А1 = (112 - 7212) 1/2 (4. 28) 

А2 = (122 - 222 )1/2 (4. 29) 


Equation (4.27) is linearized and partially difterentiated: 
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TABLE VII 


Baseline Length Deternination by Oceano 
Instruments and L. Spielvogel 


se == OG < _ му == a JG SS => == == <> «=>» cC у _к у _ и == — IUE => aa == <=> => <> => == «== AY O O O ыш щщ A < O O GEER A A O A ашк ишы == 


А 
| | 
| | 
| | 
| | 
| | Distance (п) | | 
| | | | 
| === Difference | 
| | Осеапо | Spielvogel | (п) | 
Пи ни и TAS | 
| ЕШ = Т2 | 1792.68 | 1791.7 | 0.98 { 
| go ТЗ | 1997.46 | MSDS O | 1.56 | 
| Бо = THÀ | 1973.92 | 12726 | - 1:32 
| T4 - T1 | 186 1.86 | 1860.5 | 1236 | 
| ШІ - T3 ! ОПО 2209 | 2790.8 | 1289 | 
| 12 - T4 | 2592.18 | 259073 | 1.88 | 
| EE. E 
OF 
ser -А1-11,1 | 4. 30 
ӘСІ 
OF 
BE = -A2-1 4. 
512 А L2 (4. 31) 
p = 1 (4. 32) 
95 
ӘР 
——= 41-171 (2293 
941 | | 
oF 
--- A2-172 4.34 
922 
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The Best Fit of Parabola for T1-T2 Crossing. 


Figure 4.8 
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The Best Fit of Parabola for T2-T3 Crossing. 


Figure 4.9 
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Ine Best Fit of Parabola for T3-T4 Crossing. 


Figure 4.10 
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Parabola for T1-T4Y4 Crossing. 
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Figure 4.11 
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The Best Fit of Parabola for T1-T3 Crossing. 


Figure 4.12 
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The Best Fit of Parabola for T2-T4 Crossing. 


Figure 4.13 
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Figure 4.14 Baseline Crossing Configuration. 


e 


Тае resultant observation equation for a leasr-squares solu- 


tion усл опаа 17727 р:217): 
AO + ВУ + ЫИ = 0 (4. 35) 


which can be written as: 


-/ E f O 
/ AL ZI, AZ Z2 “Al Li, - A222, O ЛЕК) 2 


AS О „О -AT Ll; Y "2 


wa пад ` |. | | 
| | о ЛЕ (4. 35) 


P LES 


Thus, the ¿east-squares solution (Mourad, 1972, p.21) gives; 
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A z-[ B, «af (Bp-1 af) 1a rna? (BP18")-18 (4. 37) 


V --p-ig* (gp-1gf) - 1 (AA +i) (4.38) 
The variance of unit weight б, іс: 
ae t 
0, = ([(BP-1B )-! (AA*W)] И/аЕ} 1/2 (4. 39) 


The variance-covariance matrix for the adjusted parameters Xo 


d, [P, +A* (Bp-1B*)-1a]-2 (4. 4) 
The adjusted parameters Xg is: 


Хо = X +4 (4. 41) 


A is a coefficient matrix for the unknoens correction, 
B is a coefficient matrix for tne observation, 
v is a vector of residuals, representing the 
corrections to observed ranges, L1 and L2, 
^ is the correction to assumed X,, 
W is a constant matrix in Equation (4.27), when 
using assumed value Xy (So «L1,¿+L2o) e 
P, is the weignt matrix of L,, 
P, is the weighting function associated with X,, 
6, is the variance of unit weight 
Xa 1S true value of 5, 21 апа 22, апа 


X, 1S an approximate value of Xg- 


The weight matrices are taken to be unity as the 


measurements were of equal precision. Table VIII shows the 
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resulting horizontal distances between transponders and the 
depths of each transponder. These results are obtained by 
measuring pairs of slant ranges at the instant the ship 
crosses the baseline. Since the ship crossed the baseline 
four times for the baseline crossing method, there are only 
four pairs of slant ranges which can ke used; these repre- 
sent the smallest data sets needed for a least squares solu- 
tion. 

The baseline lengths computed by Oceano Instruments 
and Dr. Spielvogel are given in Table IX The difference 
between Oceano's and Spielvogel's computed values are within 
the precision of the system. My results in Tabie VIII 
differ from TableIX,since Oceano and Spielvogel used much 
larger data sets (more than 250 data points), whereas I used 
the much smaller sets (11 points for each of four crossings) 
required by the two baseline crossing uethods, namely eleven 
data points for each of řour crossings. Since Oceano and 
Spielvogel used a resection method with four transponders, 
it is not necessary for tae ship to travel the baseline 
repeatediy: The only requirement is that the ship be in the 
vicinity of the center of the experimental area, where data 
can be collected continuously. Оп the contrary, the 
crossing methods require a large amount of ship time to 
collect a small amount of data. Thus, it is not surprising 
that there should be good agreement between Oceano and 
Spielvogel and rather poorer agreement between my resuits 


and theirs. 
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| | 
| TABLE VIII | 
| Results of Modified Baseline Crossing Method | 
mr | 
| Trans- Horizontal | Transponder Mean | 
| ponder Distance Between { Number Depths (п) | 
| Pairs Transponders (m)| j| 
! 1-2 1782.9119.51 | 1 1361, 5244.53 | 
| 2-3 1993-0 2560.20 | 2 133999) 0]*2 258 l 
| 3-4 1964. 34+8. 32 | B 1375: 70:5. 21 | 
| 1-4 1836-6549. 21 | 4 1407. 08+3.47 | 
| 4-3 2815.47+15.60 | | 
| 2-4 25714. 13:7 555 | | 
AE 1] D | 


6 1 


Do prueva Emo as lero c M Lum o e o lio PS Ak oe Ge OS бањи, Gmo aaa пио ИЙ 


Lines Oceano 
Т1-Т2 1792.68 
Т2-Т3 1997.46 
13-74 1913.92 
T4-T1 1861.86 
о 
T2-T4 2590 


i" ipen» Er pen been TER fo TE hjå a M a’ M a E aiaa Enn Gr EE) jigg EA Gir EE hj == EE EE GE, ED 


TABLE IX 


Spielvogel 


TUR 
193529 
1972.6 
1860.5 
2190.8 
2590.3 
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Comparison of Results for Different 


Basic 
Baseline 


Crossing 


1792.84+1. 14 
1992.42 40935 
1970-54+1. 12 
1859.00+1.04 
2789.76%0.89 
2588.74+0.82 


Methods 


Modified 
Baseline 


Crossing 


1782.91%9.51 
199 320 246250 
196 4.3448. 32 
1836-65+9.21 
2815.47+15.6 
2574.13+7.55 


| „ењ. su faro ana Mon 


V. CONCLUSIONS AND RECOMMENDATIONS 

Two methods were used to determine the horizontal 
distances between the elements of an array of four acoustic 
transponders lying on the ocean bottom at an average depth 
of 1370 meters, namely basic and modified baseline crossing 
methods. The principal differance between the two tech- 
niques is that the former requires precise depth informa- 
tion, whereas the latter does not. A "cloverleaf" maneuver 
was described which may be used to determine the depths of 
Bilemirdansponders. Because a larger number of acoustic range 
observations was used for the basic method, the standard 
deviations for the horizontal distances between transponders 
using that technique were about one-tenth those found using 
the modified method. Our results for the basic method are 
consistent with independent calculations made by Spielvogel, 
who used a much larger data set. The differences between 
our results using the modified method and Spielvogel's 
are within our relatively large standard deviations, with the 
exception of the distance measured between transponders I 
and 3. This discrepancy is probably due to the very oblique 
crossings of the 1-3 baseline which were used. 

In addition to the smallness of the data sets used for 
our calculations other sources of error in our results may 
stem from the variabilitv of the ship's speed and its pitch 
Ро ||, which we did*not take into account. 

If the absolute positions of the bottom transponders are 
required, the position of the ship may be determined by an 
absolute positioning system such as GPS while measuring 
Slant ranges from the ship to each transponder. When at 
least three or more absolute positions of the ship's 


acoustic transducer are known, the transponder coordinates 
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may then be adjusted by the method of least sguares to 
obtain a best fit between measured and calculated slant 


ranges from the known transducer positions. This method 
introduced in the Appendix. 


15 
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APPENDIX A 
ABSOLUTE POSITION OF OCEAN BOTTOM TRANSPONDERS 


The absolute position of a ship can be obtained from the 
Global Positioning System (GPS). Here, ship position деап5 
the position of tne GPS antenna. The antenna position is a 
function of time. Once the absolute position of the antenna 
is obtained, and Simultaneously, the ship's speed, heading, 
pitch, roll and acoustic data, then the absolute position of 
antenna can be used to calculate the absolute position o£ 
the bcttom transponders. The procedures are: | 

1. To deternine the absolute position of the GPS 

antenna; 

2. To convert the absolute position of antenna to the 
hull-mounted transducer using the ship heading, pitch 
and roll data; and 

3. simultaneously, to convert the Ship transducer's 
position to the pottom transponders by using acoustic 
data, 

A method 15 introduced here to deternine the absolute 
position of sea-floor transponders, if the absolute position 
of surface transducer is given.  Least-squares adjustments 
must also be supplied to get the best fit or the transponaer 
position. 

The three-dimensional geometry used is iilustrated in 
FIgure A.l. There are 4 transponders (T1, Т2, T3 and T4). 
Their estimated coordinates are (X:,Y; ,2;), the absolute 
position of ship-mounted transducer is 5(0;,7:.4,), апа «пе 
four slant ranges are 2.. 

The adjusted distances l, between transducer anrd tran- 
Sponders can be computed given coordinates of any two 
points: 
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k3? 


=> 


қ; 


T, 


[ (0-Х1) 2+ (У-Ү1) 2+ (9-21) 2 ]1/2 (А. 1) 
[ (0-Х2) 2+ (У-Ү2) 2+ (9-22) 2 41/2 (А2 2) 
[ (U-X3) 2* (V-13)2* (W-23) 2 ] 17? (A.3) 
[ (U-X4) 2+ (V-14) 2* (6-24) 2 ]1/2 (A. 4) 


which can be written simply as follows: 


7, = [ (0j WE (У: -Ү; }.2 + (1. BOE 1/2 (A.5y 


where 1 


Us V, $ 
REY, 7 
7. 


L 


are 1,2,3,4. index corresponding to 
each transponder; 

are 1,2 өе М. transducer positions 
at the surface; 

are coordinate of transducer; 

are coordinate of transponder; and 


are adjusted distance of slant range. 


It was necessary to determine the minimum number of 


positions (N) ot the null-mounted transducer that were 


needed to fix the position of four transponders in the 


taree-dimensional coordinate system. There are three 


unknown variabies to be solved for each transducer position, 


and there are four transducer positions to solve for. To 


Solve for tnese twelve unknown variables and apply ieast- 


Squares adjustment, taie number of equations must be greater 


than cr equal to the number of unknown variables. Tais is 


given by: 
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UN > 12 (А.6) 


where Ч is four equations for each interrogation; 
N is number of interrogation; and 


12 is twelve unknown for 4 transporders coordinate. 


The solution of Equation (A.6) is N23. That is, the 
absolute positions in three-dimensional coordinates of all 
transfonders can be determined with as few as three consecu- 
tive range measurements from four fixed transponders. 

Mikhail states the two-dimensional distance condition 
and its linearization (Mikhail, 1981, pp266-268), which can 
be extented to three-dinensionai coordinates. - Equation 


(A. 5) can pe linearized: 
1=10+[ (520) /19 1ьх+Е (521) 20 Аз (2-3) 749 4 ° А.Т) 
The adjusted distance is then: 
7 o (1.8) 


чһеге 2. are observed valuesof the distance; and 


у, аге corresponding residuals. 


Equation (A.7) then becomes: 


v-[ (X-U)/29 JaX-[ (Y-V) /19 Jav-[ (273) /29 ]az- 29- Q. (2.9) 
L L E ( L L ( C С C L 
Letting 
Е = до - 0. (А. 10) 


in matrix form Equation (A.9) becomes: 
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Vi Ба 6,2 bj3 O 0 x AX, E 
Vz 0 0 0 Dai b22 b23 O > aY, i 
V3 > е . • å ` з АЯ, К s 
t d бе 11 
Y % o o % е | | B ta 
| ЕЕ ЛАШ | (М. ali 
VEE 
which may be concisely written as: 
a + BW x e = | (A. 12) 
where 


V is residual of the observation; 
B is numerical coefficient of tne unknown 
variable correction ; 
AX is unknown variable correction ; апа 


f is numerical constant term. 


The weight matrix can be derived from the standard devi- 
ation of the slant range, which discussed in the previous 
ШОО рег. The least squares solution of Equation (A.12) is 
(Mikhail, 1981, p259): 


NES BC WB (А. 13) 
t = B WE (А. 14) 
ELI uot (A. 15) 
X -X, tax (A. 16) 
v = f - ВАХ (A. 17) 


o9 


oD 


N is 
Я 15 
t is 
Ах 15 
Xg ls 
^ . 
X is 
б; LS 
„~ 
02 
о 
Г 15 


15 


N is 


3 15 


-[ G2 (B N-1 ве) ]1/2 


V NW v / T 


(A. 19) 


{А ЖЭ) 


(А. 20) 


coefficient matrix of the normal equations, 


weight matrix, 


vector 
vector 
vector 


vector 


of 'constants' in the normal equations, 
Of parameter corrections, 
of approximate vaiue, 


of adjustea value, 


standard deviation of adjusted distance, 


estimate oí the reference variance, 


number 
nunber 


number 


problem. 


of statical degrees of freedom, 
of given transducer position, and 


of observations necessary to specify 


uniquely tne nodel that underlies the adjustuent 


Since during linearization we neglect all second and 


higher order terns, 


we must ensure that X is not signifi- 


cantly іп error because of this approximation. 


Conseguently, 


Rew approximation and compute anotaer correction. 


we must iterate the solution by using X ās a 


This 


procedure is repeated until the correction is insignificant. 


rhe final result is: 


where 


(A. 21) 
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n is total iteration, 
а is adjusted value at n iteration, 
EC is correction at n iteration can be neglected compare 
with previous correction or within tolerance 
Correction, and 


Ха 15 true value of each variable. 


Since errors in transducer's position exist, it now 
must te established how this error affects tha accuracy of 
transponder positions. This is discussed below: 

Let the standard deviaticns of transducer coordinates be 
+0 y ТО +0 - Кеса11 Eguation (A.5): 


о 


7. = [(U -X, )2+ (V -Y, ) 2*(N 72; ) 2 j172 TOD 


The unknown coordinates.X,Y,Z , are: 


= Uus [ 22- (V-Y.) SL JAS Е 

у. = У - [32- (9-2. ) 2- (0-Х. ) 2] 1/2 (A. 24) 

L L L » 

л: Д [ Y2- (U- X, ) 2- (V- Y, ) 2317? (А. 25) 
C L 


Now assume 


Е1 = [ 22- (V-Y; ) 2- (W-2;)2]?/2 (A. 26) 

F2 = [ 22-(W-2.) 2-(U-X.) 2]2/2 (A. 27) 
t L . с 

ЕЗ - Г22-(0-Х,) 2-(У-Ү.)211/2 (A. 28) 
E L e 


So, the variances of Equations (A.23), (A.24) and (A-25) are 
(Mikhail, 1981, p.181): 


дг =02-(0; /Е1) 2024 (V-Y, /F1)202+ (8-2; /F1) 202 (A. 29) 
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de =02-(1. /Р2) 202+ (5-7, /Р2) 202+ (0-Х, 7/72) 202 (A. 30) 


02 =02- (1. /ЕЗ) 242+ (0-Х, /Р3) 202+ {У-1, /РЗ) 202 (4.31) 
and the position accuracy is: 


0, = + (62 +02 БОЈЕ ) 1/2 (А. 32) 
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